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Abstract 

We analyze the Landau-Lifshitz-Gilbert equation when the precession motion of the magnetic 
moments is additionally subjected to an uniaxial anisotropy and is driven by a multiplicative cou- 
pled stochastic field with a finite correlation time r. The mean value for the spin wave components 
offers that the spin-wave dispersion relation and its damping is strongly influenced by the deter- 
ministic Gilbert damping parameter a, the strength of the stochastic forces D and its temporal 
range r. The spin-spin-correlation function can be calculated in the low correlation time limit by 
deriving an evolution equation for the joint probability function. The stability analysis enables us 
to find the phase diagram within the a — D plane for different values of r where damped spin wave 
solutions are stable. Even for zero deterministic Gilbert damping the magnons offer a finite life- 
time. We detect a parameter range where the deterministic and the stochastic damping mechanism 
are able to compensate each other leading to undamped spin-waves. The onset is characterized by 
a critical value of the correlation time. An enhancement of r leads to an increase of the oscillations 
of the correlation function. 
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I. INTRODUCTION 



Magnetism can be generally characterized and analyzed on different length and time scales. 
The description of fluctuations of the magnetization, the occurrence of damped spin waves 
and the influence of additional stochastic forces are successfully performed on a mesoscopic 
scale where the spin variables are represented by a continuous spatio-temporal variable [l]. 



In this case a well established approach is based upon the Landau-Lifshitz equation [2j which 
describes the precession motion of the magnetization in an effective magnetic field. This 
field consists of a superposition of an external field and internal fields, produced by the in- 
teracting magnetic moments. The latter one is strongly influenced by the isotropic exchange 
interaction and the magnetocrystalline anisotropy, for a recent review see The studies 
using this frame are concentrated on different dynamical aspects as the switching behav- 
ior of magnetic nanoparticles which can be controlled by external time-dependent magnetic 
fields j^j and spin-polarized electric currents {5, (J. Such a current- induced spin transfer 
allows the manipulation of magnetic nanodevices. Recently, it has been demonstrated that 
an electric current, flowing through a magnetic bilayer, can induce a coupling between the 
layers [7|. Likewise, such a current can also cause the motion of magnetic domain walls in 
a nanowire [S|. Another aspect is the dynamical response of ferromagnetic nanoparticles 
as probed by ferromagnetic resonance, studied in Q]. In describing all this more complex 
behavior of magnetic systems, the Landau-Lifshitz equation has to be extended by the in- 
clusion of dissipative processes. A damping term is introduced phenomenologically in such a 
manner, that the magnitude of the magnetization S is preserved at any time. Furthermore, 
the magnetization should align with the effective field in the long time limit. A realization 
is given by Q 

— = - 7 [SxB cff ]-£[Sx (SxB cff )] . (1) 
The quantities 7 and e are the gyromagnetic ratio and the damping parameter, respectively 



An alternative equation for the magnetization dynamics had been proposed by Gilbert [10]. 
The Gilbert equation yields an implicit form of the evolution of the magnetization. A com- 
bination of both equations, called Landau-Lifshitz-Gilbert equation (LLG) will be used as 
the basic relation for our studies, see Eq. (J2J). The origin of the damping term as a non- 



relativistic expansion of the Dirac equation has been discussed in ll| and a generalization 



of the LLG for conducting ferromagnetics is offered in 



12j . The form of the damping seems 



to be quite general as it has been demonstrated in 13| using symmetry arguments for fer- 
roelectric systems. 

As a new aspect let us focus on the influence of stochastic fields. The interplay between 



current and magnetic fluctuations and dissipation has been studied recently in [14|. Via the 
spin-transfer torque, spin-current noise causes a significant enhancement of the magnetiza- 
tion fluctuations. Such a spin polarized current may transfer momentum to a magnet which 
leads to a spin-torque phenomenon. The shot noise associated with the current gives rise to 
a stochastic force [15(. In our paper we discuss the interplay between different dissipation 
mechanism, namely the inherent deterministic damping in Eq. ([I]) and the stochastic mag- 
netic field originated for instance by defect configurations giving rise to a different coupling 
strength between the magnetic moments. Assuming further, that the stochastic magnetic 
field is characterized by a finite correlation time, the system offers memory effects which 
might lead to a decoherent spin precession. To that aim we analyze a ferromagnet in the 
classical limit, i.e., the magnetic order is referred to single magnetic atoms which occupy 
equivalent crystal positions, and the mean values of their spins exhibit a parallel orientation. 
The last one is caused by the isotropic exchange interaction which will be here supplemented 
by a magneto-crystalline anisotropy that defines the direction of the preferred orientation. 
Especially, we discuss the influence of an uniaxial anisotropy. The coupling between differ- 
ent dissipation mechanisms, mentioned above, leads to pronounced correlations, which are 
discussed below. Due to the multiplicative coupling of the stochastic field and the finite 
correlation time the calculation of the spin-spin correlation function is more complicated. 
To that aim we have to derive an equivalent evolution equation for the joint probability 
distribution function. Within the small correlation time limit this app roach can be fulfilled 
in an analytical manner. Our analysis is related to a recent paper [16| in which likewise the 
stochastic dynamics of the magnetization in ferromagnetic nanoparticles has been studied. 
Further, we refer also to a recent paper [l7| where the mean first passage time and the 
relaxation of magnetic moments has been analyzed. Different to those papers our approach 
is concentrated on the correlation effects in stochastic system with colored noise. 
Our paper is organized as follows: In Sec. II we discuss the LLG and characterize the ad- 
ditional stochastic field. The equations for the single and the two particle joint probability 
distribution are derived in Sec. III. Using these functions we obtain the mean value of the 
spin wave variable and the spin-spin correlation function. The phase diagram, based on the 
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stability analysis, is presented in Sec. IV. In Sec.V we finish with some conclusions. 



II. MODEL 

In order to develop a stochastic model for the spin dynamics in ferromagnetic systems let 
us first consider the deterministic part of the equation of motion. We focus on a description 
based upon the level of Landau-Lifshitz phenomenology [2|, for a recent review see [3]. To 
follow this line we consider a high spin systems in a ferromagnet sufficiently below the 
Curie temperature. In that regime the dynamics of the magnet are dominated by transverse 
fluctuations of the spatio-temporal varying local magnetization. The weak excitations, called 
spin waves or magnons, are determined by a dispersion relation, the wavelength of which 
should be large compared to the lattice constant a, i.e., the relation q ■ a <C 1 is presumed to 
be satisfied, where q is the wavenumber. In this limit the direction of the spin varies slowly 
while its magnitude |S| = m s remains constant in time. A proper description for such a 



situation is achieved by applying the Landau-Lifshitz-Gilbert equation (LLG) [4j, [10|, Il8 |. 
The spin variable is represented by S = m s n, where n(r, t) is a continuous variable which 
characterizes the local orientation of the magnetic moment. The evolution equation for that 
local orientation reads 

^ = -^nx[B cS + a[iixB cf ]]. (2) 

The quantities 7 and a are the gyromagnetic ratio and the dimensionless Gilbert damping 
parameter, respectively, where a is related to e introduced in Eq. flTJ. B c g is the effective 
magnetic field that drives the motion of the spin density. Generally, it consists of an internal 
part originated by the interaction of the spins and an external field. This effective field is 
related to the Hamiltonian of the system by functional variation with respect to n 

Bo^-^f • (3) 
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In absence of an external field the Hamiltonian can be expressed as 

H = d 3 r {w ex + w an } , with 

11 (4) 

w ex = -m s K(Vn) 2 and w an = ^m s T sm 2 9 . 

Thereby, the constants k and T denote the exchange energy density and the magneto- 
crystalline anisotropy energy density. To be more precise, k oc Ja 2 , J being the coupling 



strength that measures the interaction between nearest neighbors in the isotropic Heisenberg 



model 21] . Once again a is the lattice constant. Notice that the form of the exchange energy 
in the Hamiltonian (j3J) arises from the Heisenberg model in the classical limit. The quantity 
9 represents the angle between n and the anisotropy axis i> = (0,0,1), where v points 
in the direction of the easy axis in the ground state in the case of zero applied external 
field. Thus, the constant T > characterizes anisotropy as a consequence of relativistic 
interactions (spin-orbital and dipole-dipole ones |20j). In deriving Eq. (j4]) we have used 
n 2 = 1. Although it is more conventional to introduce the angular coordinates (6, $) 

aa, 

we find it more appropriate to use Cartesian coordinates. To proceed, we divide the vector 
n into a static and a dynamic part designated by [i and (p, respectively. In the linearized 
spin wave approach let us make the ansatz 

n(r,t) = fj,(r) + (p(r,t) — fj,i> + <p , = const. , (5) 

where n 2 = 1 is still valid. The effective field can now be obtained from Eqs. fl3]) and (j4j). 
This yields 

B eff = KV>-rV; <S = {<p 1 ,<p 2 ,0). (6) 

Eq. (j2j) together with Eqs. ([3]) and (jlj) represent the deterministic model for a classical 
ferromagnet. In order to extent the model let us supplement the effective magnetic field in 
Eq. (jSJ) by a stochastic component yielding an effective random field B eff = B eff + rj(t). The 
stochastic process rj(t) is assumed to be Gaussian distributed with zero mean and obeying 
a colored correlation function 



t-f 



T i.i 



(7) 



Here, Dij and are the noise strength and the finite correlation time of the noise r\. 
Due to the coupling of the effective field to the spin orientation n the stochastic process 
is a multiplicative one. Microscopically, such a random process might be originated by 
a fluctuating coupling strength for instance. The situation associated with our model is 
illustrated in Fig. [TJ and can be understood as follows: The stochastic vector field ry(t) is able 
to change the orientation of the localized moment at different times. Therefore, fixed phase 
relations between adjacent spins might be destroyed. Moreover, the rj(tk) are interrelated 
due to the finite correlation time r. The anisotropy axis defines the preferred orientation of 
the mean value of magnetization. Due to the inclusion of r](t) the deterministic Eq. (J2J) is 



anisotropy axis 



random field at 
different times U 




exchange oc J 



FIG. 1. Part of a ferromagnetic domain inffuenced by stochastic forces for the example of cubic 
symmetry with lattice constant a. The black spin in the center only interacts with its nearest 
neighbors (green), where J is a measure for the exchange integral. 



transformed into the stochastic LLC Using Eq. 

dcp 7 
~dt = 



it follows 

H + (p)x [B eff + a [O + <p) x Bcff]] . 



1 + a 2 

The random magnetic field is defined by 

k V 2 (f - T <p' + rj(t) , 



B 



cff 



(8) 



(9) 



where <p' is given in Eq. ([6]). With regard to the following procedure we suppose the random 
field to be solely generated dynamically, i.e., n x rj(t) = ip x f](t). So far, the dynamics 
of our model (Eqs. © and ([9])) are reflected by a nonlinear, stochastic partial differential 
equation (PDE). Using Fourier transformation, i.e., ip(c\, t) — T{<p(r, t)} and introducing 
the following dimensionless quantities 



K 
f 
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7 r 



t — ut 



X(t) 



the components ^j(q, t) fulfill the equation 

d 



dt 



(10) 



(11) 



The quantity Zo is the characteristic magnetic length 22]. The vector Q and the matrix A 
are given by 



'-(a/^i + V^ 

ipi - a/j, ip 2 




1 + a 2 



(12) 



and 



-ip 3 afj,ip 3 ipt- afj,ip 2 

1p 2 -Ipl 



(13) 



For convenience we have substituted t — > t again. The statistical properties of \(t) are 
expressed as (A(t)) = and 



D 



Xki(t,t') = (X k (t) Xtf)) = —5 kl exp 

Tkl 



t-f 

Tkl 



¥2D kl 6 kl 6(t-t') 



(14) 



Incidentally, in the limit r — > the usual white noise properties are recovered. We empha- 
size that although we regard the long- wavelength limit (a • q 1), wave vectors for which 
lo ■ q 1 (in Eq. ( FlOl) ) can also occur [22|. But this case is not discussed in the present 
paper and will be the content of future work. Whereas, in what follows we restrict our 
considerations to the case q — > so that, actually, Iq ■ q 1 is fulfilled. Hence, we can set 
(3 = 1 approximately in Eq. (flQj) . Due to the anisotropy the spin wave dispersion relation 
offers a gap at q = 0. Owing to this fact if) is studied at zero wave vector. For this situation 
the assumption of a space-independent stochastic force ^(t), compare Eq. ([7]), is reasonable. 
For non-zero wave vector the noise field should be a spatiotemporal field r]i((r,t). Because 
our model is based on a short range interaction we expect that the corresponding noise 
correlation function is 5-correlated, i.e. instead of f lT4|) we have 



Xki(r,t; r',t') = —5 M exp 

Tkl 



t-f 

Tkl 



2M5{r 



where M is the strength of the spatial correlation. Using this relation we are able to study 
also the case of small q which satisfies Iq ■ q <g 1. In the present paper we concentrate on 
the case of zero wave vector q = 0. 



III. CORRELATION FUNCTIONS 



In the present section let us discuss the statistical behavior of the basic Eqs. (lli p -( J14p . They 
describe a non-stationary, non-Markovian process attributed to the finite correlation time. 
Due to their common origin both characteristics can not be analyzed separately. In the 
limit t — > 0, Eq. ( II ip defines a Markovian process which provides also stationarity by an 



appropriate choice of initial conditions [231 ] . However, the present study is focused on the 
effect of nonzero correlation times. To that purpose we need a proper probability distribution 
function which reflects the stochastic process defined by Eqs. (jlip -( )14p . In deriving the 
relevant joint probability distribution function we follow the line given in [24|. where the 
detailed calculations had been carried out, see also the references cited therein. In particular, 
it has been underlined in those papers that in order to calculate correlation functions of type 
(ijji(t) ipj(t')) a single probability distribution function P(ifj,t) is not sufficient. Instead of 
that one needs a joint probability distribution of the form P(if), t; if)', t'). Before proceeding 
let us shortly summarize the main steps to get the joint probability distribution function. 
To simplify the calculation we assume tm = r Ski and P>ki = D Ski- Notice that our system 
has no ergodic properties what would directly allow us to relate the stochastic interferences 
with temperature fluctuations by means of a fluctuation-dissipation theorem. Based on 
Eq. (jTTj) the appropriate joint probability distribution is defined by [24j, |25|, for a more 



general discussion compare also 26]: 



P{4>, t; V', = (5(if>(t) - V) <W) - </>')> • (15) 



Here the average is performed over all realizations of the stochastic process. In defining the 
joint probability distribution function we follow the convention to indicate the stochastic 
process by the function if){t) whereas the quantity without arguments if) stands for the 
special values of the stochastic variable. These values are even relalized with the probaility 
P(ift,t]if)',t'). The equation of motion for this probability distribution reads according to 
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24] 



d_ 

Ot 



d 




Xjk(t, t\ 



XjkV-i h 



dip', 
o 

where Novikov's theorem 



(16) 



27] has been applied. Expressions for the response functions 
5ipi(t)/5 Afc(ti) and 5 ipi(t') j '5 \ k (ti) can be found by formal integration of Eq. ( ITT]) and 
iterating the formal solution. After a tedious but straightforward calculation including the 
computation of the response functions to lowest order in (t — ti) and (f — ti) and the 
evaluation of several correlation integrals referring to xu from Eq. (fT4"|) . Eq. ( fl6|) can be 



rewritten in the limit of small correlation time r as 
^P s (i>,t;i>' 1 t>)={C°(i> 1 T) 



+ exp [_( t _ t ')/r] D — A lk {xp) -^r KkW) | Ps(fl>, t; </>', t') . 



(17) 



Thereby, transient terms and terms of the form oc r exp[— (t — t')/r] (these terms would lead 
to terms of order r 2 in Eq. ( 122]) ) have been neglected. The result is valid in the stationary 



case characterized by t — > oo and if — Y oo but finite s = t — t'. In Eq. ([T 
appearing in the equation for the single probability density. Following 2 
reads 



£° is the operator 



the operator 



d d d 

fyi lKV> lkK ^'dip n 



D [A nk (^)-rM nk (^)] 



+ D 2 t 



d Id 
K nkm (ip>)—A lm (xp>) + -A nrn {il})—Ki km {tl)) 
oipi 2 d% 



(18) 



with 



M 



nk 



dA nk _ dtt r 



TlZfc 



A 



dlp r 

dA nl 



dip. 
9A nk 



(19) 



A 



rl ■ 



dlp r dlp r 

The equation of motion for the expectation value (ipi) s can be evaluated from the single 
probability distribution in the stationary state 

d 



at 



P s (iP,t) = £ P s (iP,t). 



(20) 



One finds 
d 



<)X,,k (A nk -rM nk ) 



D 2 



d fdA 



dip r \ dip. 



^ik 



^nkra I A r7y 



d (dA lk 



*" 2 \ cfyV V <^n 



A., 



nm I -L*-rkm 



(21) 



The knowledge of the evolution equation of the joint probability distribution P(ip,t; ip',t') 
due to Eqs. f lT7|) and ffT8l) allows us to get the corresponding equation for the correlation 
functions. Following again [241], it results 



ik 



dipr, 



(A nk - r M nfc ) 



D 2 r 



d fdki 



ik 



dlp r \dlpr, 

d ( dA ik 



nkm I !*-rm 



A, 



.D exp 



dip r \ dip, 

t-t' 



-Ar. 



nm I K r km 



(22) 



T 



(A ik (^(t))A jk (^(t'))) t 



where the symbol [...] t denotes the quantity [...] at time t. As mentioned above the result 
is valid for t, t' — > oo while s = t — t' > remains finite. The quantities M nk and K k i m are 
defined in Eq. ( Tl9l) . The components fli and A;.,- are given in Eqs. f fl~2l and ( TTBl . Performing 
the summation over double-indices according to Eqs. (12Tj) and (122|) we obtain the evolution 
equations for the mean value and the correlation function 

d 



dt 



(1>i(t)). = G* (Mt)) s i 



(23) 



and 



Ts c « (s) 



ds 



(^.(f + 8 ) ip j ( t ')) s =G lk (ip k (f + s) iP,(t')) f 



+ D exp 



T 



(A ik (i/>(t' + s))A jk (W))) s . 



(24) 



Notice, that in the steady state one gets Cij(t,t') = Cij(s) with s = t — t'. The matrix 
components of Gi k are given by 



Gi k 



-Ay A 2 
-A 2 -Ay 

-A 3 



(25) 
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where 

A x = -D 2 r(6fi 2 a 2 - 1) £ 4 + 2fx 2 aDr £ 3 - D(/i 2 a 2 - 2) £ 2 + fi 2 a £ 

A 2 = ^ftaD 2 r (11 - 3/iV) £ 4 + hDt (/iV - l) £ 3 + 3fiDa£ 2 - /x£ (26) 

A 3 = +D 2 r (3fi 2 a 2 + 1) £ 4 - V«Dr^ + 2D £ 2 , 

and £ is denned in Eq. (|T2|) . At this point let us stress that in the case £' = the term 
oc exp[— (£ — £') /r] on the rhs. in Eqs. ( |22|) and ( l24l) . respectively, would vanish in the steady 
state, i.e. 

(Mt' + s)^(t')) s ^(^(s)^(0)) s . 

The occurrence of such a term is a strong indication for the non-stationarity of our model. An 
explicit calculation shows, that in general this inequality holds for non-stationary processes 
23]. 



IV. RESULTS 



The solution of Eq. (1231) can be found by standard Greens function methods and Laplace 
transformation. As the result we find 

e~ Alt cos(A 2 £) e~ Al 1 sin(A 2 £) ^ 
-e~ Alt sin(A 2 £) e~ Al 1 cos(A 2 £) 



V 



o 



o 



-A 3 t 



(V>o> f 



(27) 



where (V>o) s = (VK^ = 0)) s are the initial conditions. The parameters Ax, A 3 and A 2 defined 
in Eqs. ( )26l) play the roles of the magnon lifetime and the frequency of the spin wave at 
zero wave vector, respectively. As can be seen in Eq. ( |26i) all of these three parameters are 
affected by the correlation time r and the strength D of the random force. Moreover, the 
Gilbert damping parameter a influences the system as well. The solution of Eq. (I24I) for 
the correlation function in case of £' = is formal identical to that of Eq. (|2"Tj) . The more 
general situation £' ^ allows no simple analytic solution and hence the behavior of the 
correlation function C(s) is studied numerically. In order to analyze the mean values and 
the correlation function let us first examine the parameter range where physical accessible 
solutions exist. In the following we assume (^i(O)) = (^2(0)) = (V>o) an d (^(O)) = 0, since 
the solutions for ipx(t) and ip2(t) on the one hand and ips(t) on the other hand are decoupled 
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in Eq. (|27|) . Therefore, spin wave solutions only exists for non-zero averages (ipi(t)) and 
("02 (t))- The existence of such non-trivial solutions are determined in dependence on the 
noise parameters D and r and the deterministic damping parameter a. Notice, that the 
dimensionless quantity D = D/T, i.e., D is the ratio between the strength of the correlation 
function (Eq. ([7])) and the anisotropy field in the original units. The stability of spin wave 
solutions is guaranteed for positive parameters A\ and A3. According to Eqs. (1261) the 
phase diagrams are depicted in Fig. [2] within the a — D plane for different values of the 
correlation time r. The separatrix between stable and unstable regions is determined by 
the condition Ai = 0. The second condition A 3 = is irrelevant due to the imposed initial 
conditions. As the result of the stability analysis the phase space diagram is subdivided into 



four regions where region IV does not exist in case of r = 0, see Fig. 2(a) For generality, 
we take into account both positive and negative values of D indicating correlations and 
anti-correlations of the stochastic field. Damped spin waves are observed in the areas I and 
IV, whereas the sectors II and III reveal non-accessible solutions. In those regions the spin 
wave amplitude, proportional to exp[— Ait], tends to infinity which should not be realized, 



compare Figs. 2(b) 2(d) Actually, a reasonable behavior is observed in regions I and IV. As 
visible from Fig. [2] damped spin waves will always emerge for D > even in the limit of zero 
damping parameter a and vanishing correlation time r. This behavior is shown in Fig. [3l 



where the evolution of (?/>i(t)) is depicted for different values of a. As can be seen in Fig. 2(a" 



the solution for D < is unlimited and consequently, it should be excluded further. Contrary 
to this situation, additional solutions will be developed in region IV in case of r > and 



simultaneously a = 0, see Figs. 2(b) 2(d) Thereby the size of area IV grows with increasing 
t. Likewise, the extent of region I decreases for an enhanced r. However, in the limit of 
D = and consequently for r = 0, too, only damped spin waves are observed. Immediately 
on the separations line undamped periodic solutions will evolve, compare the sub-figures in 
Fig. [2l This remarkable effect can be traced back to the interplay between the deterministic 
damping and the stochastic forces. Both damping mechanism are compensated mutually 
which reminds of a kind of resonance phenomenon. The difference to conventional resonance 
behavior consists of the compensation of the inherent deterministic Gilbert damping and the 
stochastic one originated from the random field. This statement is emphasized by the fact 
that undamped periodic solutions do not develop in the absence of stochastic interferences, 
i.e., D = 0. The situation might be interpreted physically as follows: the required energy 
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FIG. 2. a — D plane for fixed magnetization [i = 0.9 and different values of r. 



that enables the system to sustain the deterministic damping mechanisms is delivered by 
the stochastic influences due to the interaction with the environment. To be more precise, in 
general, the Gilbert damping enforces the coherent alignment of the spin density along the 
precession axis. Contrary, the random field supports the dephasing of the orientation of the 
classical spins. Surprisingly, the model predicts the existence of a critical value r = r c > 
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1.00-1 




FIG. 3. Evolution of the mean value (ijji(t)), with = 0.9, D = 0.1 and r = 0. a varies from 
(dash-dotted line), 0.05 (solid line), 0.5 (dotted line) and 1 (dashed line). 

depending on a and D which determines the onset of undamped periodic solutions. Notice, 
that negative values of r c are excluded. The critical value is 

[fi 2 (a 3 -Da 2 + a) + 2D](l + a 2 f 
Tc ~ 2D ^ (a 3 - 3Da 2 + a) + D 2 ' [ ' 

Hence, this result could imply the possibility of the cancellation of both damping processes. 
Examples according to the damped and the periodic case are displayed in Fig. HJ An increas- 



ing t favors the damping process as it is visible in Fig. 4(a) Based on estimations obtained 



for ferromagnetic materials [29) and references therein, the Gilbert damping parameter can 
range between 0.04 < a < 0.22 in thin magnetic films, whereas the bulk value for Co takes 
ab ~ 0.005. The phase space diagram in Fig. [2] offers periodic solutions only for values of 
a larger than those known from experiments. Therefore such periodic solutions seem to be 
hard to see experimentally. We proceed further by analyzing the behavior of the correlation 
function by numerical computation of the solution of Eq. (|24p with Eqs. 0251) and (1261) . As 
initial values we choose Cik(t = t', t') = Cik(s = 0) = Cq for every combination i,k = {1, 2, 3}. 



The results are depicted in Figs. [5] and El Inspecting Figs. 5(a) 5(c) one recognizes that an 



enhancement of the correlation time r leads to an increase of the oscillations within the 
correlation functions Cik, k = {1,2,3}. Moreover, Fig. 5(d) reveals that the oscillatory 
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t t 

(a) (b) 

FIG. 4. Evolution of the mean values (ipipit)), with /i = 0.9. (a): D = 0.1, a = 0.005 and r varies 
from 10 (solid line), 1 (dotted line) and (dash-dotted line), (b): D = 2, a = 1 and r = r c 1.79 
(Eq. ([28j) b The solid line represents and the dash-dotted line is (V^)- 



behavior of C31 seems to be suppressed. Obviously, the decay of the correlation function is 



enhanced if r growths up. The pure periodic case for r = r c , corresponding to Fig. 4(b) 



is depicted in Fig. |6j Exemplary, C 12 and C31 are illustrated. The behavior of the latter is 



similar to the damped case, displayed in Fig. 5(d) , unless slight oscillations occur. However, 



if one compares the form of Cu in Fig. 5(b) and Fig. [6] the differences are obvious. The am- 
plitude of the correlation function for the undamped case grows to the fourfold magnitude 
in comparison with Co, whereas the damped correlation function approaches zero. Further, 
a periodic behavior is shown in Fig. [61 and therefore the correlation will oscillate about zero 
but never vanish for all s — t — t' > 0. 



V. CONCLUSIONS 



In this paper we have analyzed the dynamics of a classical spin model with uniaxial 
anisotropy. Aside from the deterministic damping due to the Landau-Lifshitz-Gilbert 
equation the system is subjected to an additional dissipation process by the inclusion of a 
stochastic field with colored noise. Both dissipation processes are able to compete leading to 

15 




a more complex behavior. To study this one we derive an equation for the joint probability 
distribution which allows us to find the corresponding spin-spin-correlation function. This 
program can be fulfilled analytically and numerically in the spin wave approach and the 
small correlation time limit. Based on the mean value for the spin wave component and 
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FIG. 6. Correlation functions C ik (s) for r = r c « 1.79 (Eq. flM}), H = 0.9, D = 2 and a = 1. The 
dotted line represents C12 and the solid line is C31. 



the correlation function we discuss the stability of the system in terms of the stochastic 
parameters, namely the strength of the correlated noise D and the finite correlation time 
r, as well as the deterministic Gilbert damping parameter a. The phase diagram in the 
a — D plane offers that the system develops stable and unstable spin wave solutions due to 
the interplay between the stochastic and the deterministic damping mechanism. So stable 
solutions evolve for arbitrary positive D and moderate values of the Gilbert damping a. 
Further, we find that also the finite correlation time of the stochastic field influences the 
evolution of the spin waves. In particular, the model reveals for fixed D and a a critical 
value r c which characterizes the occurrence of undamped spin waves. The different situa- 
tions are depicted in Fig. [2j Moreover, the correlation time r affects the damped spin wave 
which can be observed in regions I and IV in the phase diagram. If the parameters D and 
a changes within these regions, an increasing r leads to an enhancement of the spin wave 



damping, cf. Fig. 4(a) The influence of r on the correlation functions is similar as shown 



in Figs. 5(a) 5(c) The study could be extended by the inclusion of finite wave vectors and 



using an approach beyond the spin wave approximation. 
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